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1/-^ ' Abstract 
(N ' 

T-H Context. For applications in population synthesis, libraries of theoretical stellar spectra are often considered an alternative to template libraries 
T-H of observed spectra, because they allow a complete sampling of stellar parameters. Most attention in published theoretical spectral libraries has 
, been devoted to the visual wavelength range. 

. Aims. We present a detailed comparison of theoretical spectra in the range 1.57- 1.67pm, for spectral types from A to early M and for giants 
and dwarf stars, with observed stellar spectra at resolutions around 3000, which would be sufficient to disentangle the different groups of late 
type stars. 

O I Methods. We have selected the NeMo grids of stellar atmospheres to perform such a comparison. 

Results. We first demonstrate that after combining atomic and molecular line lists, it is possible to match observed spectral flux distributions 
with theoretical ones very well for almost the entire parameter range covered by the NeMo grids at moderate resolution in the visual range. In 
the infrared range, although the overall shape of the observed flux distributions is still matched reasonably well, the individual spectral features 
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are reproduced by the theoretical spectra only for stars earlier than mid F type. For later spectral types the differences increase and theoretical 
• • spectra of K type stars have systematically weaker line features than those found in observations. These discrepancies are traced back to stem 
_ primarily from incomplete data on neutral atomic lines, although some of them are also related to molecules. 

Conclusions. Libraries of theoretical spectra for A to early M type stars can be successfully used in the visual regions for population synthesis 
» I but their application in the infrared is restricted to early and intermediate type stars. Improving atomic data in the near infrared is a key element 

in making the construction of reliable libraries of stellar spectra in the infrared feasible. 
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1. Introduction 

Due to the rapidly increasing spectral resolution of galaxy sur- 
veys (e.g. Sloan Digital Sky Survey), modelling any galaxy 
with spectral synthesis technique requires a high spectral reso- 
lution stellar library. 

Over the last years much progress has been made in syn- 
thesis models (e.g. Pelat_1997; Leitherer et al. 1999, Moultaka 
& Pelat l^OOt Bruzu al & C hariot 2003; Le Borgne et al. 2004J 
Cid Fernandes et al. I2005> . Moderate to high spectral resolu- 
tion observations of stars in order to construct reliable tem- 
plate libraries have also been performed in the visible domain 
(e.g. STELIB of Le Borgne et al. "^OOT UVES of Bagnulo et 
al. 12003 Cou deFed of Valdes et al. 2004) and in the IR (e.g. 
Dallier et al. 119961 Meyer et al. IT998I Ivanov et al. |2^ . 
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However, the major limitation of all these libraries is the sam- 
pling of stellar parameters such as metallicity. 

One way to avoid such a difficulty is to build libraries of 
theoretical stellar spectra in order to choose the desired phys- 
ical parameters. In this sense, some extensive libraries of syn- 
thetic spectra have appeared recently for the visible range. 

Murphy & Meiksin (2004), based on Kurucz's ATLAS9 
model atmospheres (Kurucz 1993aj . have built a high resolu- 
tion (A/AA - 250000) stellar library over an extended visible 
range (3000 to 10000 A). The convection zone is treated us- 
ing the Mixing Length Theory (MLT) with the overshooting 
treatment of Kurucz (cf. Castelli et al. 1997). This library pro- 
vides spectra for 54 values of effective temperature from 5250 
to 50000 K, 1 1 values of log surface gravity from 0.0 to 5.0 and 
19 metallicities from -5.0 to 1 .0. They compared their synthetic 
library with observed spectra (STELIB library of Le Borgne et 
al. 2003 ) for the colours and the Lick indices and found in gen- 
eral a good agreement. 
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Still based on Kurucz's models, but with enhanced molec- 
ular line lists and the overshooting option switched off, Munari 
et al. J2005> also present a library of synthetic spectra, for a 
similar wavelength range (2500 to 10500 A). They use a new 
grid of ATLAS9 model atmospheres (Castelli & Kurucz 2003 1. 
The effective temperature of these spectra is contained between 
3500 and 47500 K, the log surface gravity between 0.0 and 
5.0 and the metallicity between -2.5 and 0.5. These spectra are 
computed at a resolving power of AjAA - 500000 and then 
Gaussian convolved to lower resolution (< 20 000). In con- 
trast with Murphy & Meiksin J2004t . the predicted energy level 
lines are not included into the line lists used to build this library, 
as Munari et al. ( 2005 1 favour the spectroscopic use rather than 
the photometric use of the synthetic spectra. The addition of 
these "predicted lines" permits to have a better statistical flux 
distribution, but individual wavelengths can be wrong by up to 
5%. 

We would like to point out here that usually only the lower 
lying energy levels of atoms have been determined in the labo- 
ratory, particularly for complex spectra such as those from neu- 
tral or singly ionized iron. If only those transitions were taken 
into account, the atmospheric line blanketing computed from 
such data would be severly incomplete. A lot of weak lines, 
possibly unidentified even in the solar spectrum but neverthe- 
less present, would be missed. This would lead to an overes- 
timation of the ultraviolet flux which in turn would be "com- 
pensated" by a lack of flux in the visual (Kurucz 1992 1. To 
avoid this deficiency and to improve the temperature structure 
of the model atmospheres, the spectrophotometric flux distri- 
bution, and the photometric colors requires to account for lines 
for which one or both energy levels have to be predicted from 
quantum mechanical calculations. This has been one of the 
main goals of the ATLAS9 models of Kurucz (IT9921 11993a> . 
As the theoretical predictions are accurate to only a few per- 
cent, individual wavelengths can be wrong by up to a few 100 A 
in the visual. Also the line oscillator strengths are sufficiently 
accurate merely in a statistical sense. This still permits to im- 
prove the total flux within a wavelength band of a few dozen A 
in the visual, but individual features do appear in the wrong part 
of the spectrum. For spectroscopy at higher resolutions, partic- 
ularly if the spectrum is rectified or when working within small 
wavelength bands, adding the predicted energy level lines "pol- 
lutes" the theoretical spectrum with extra "noise". This has to 
be avoided in libraries devoted to automatic fitting procedures 
or if particular line features are essential to identify a certain 
spectral type (lines predicted for the wrong wavelength will 
make this more difficult). Hence, with present atomic data, ei- 
ther choice is only a compromise solution. 

By combining three different model atmospheres, the high- 
resolution stellar library of Martins et al. ( 2005 1 provides the 
largest coverage in effective temperature (from 3000 to 50000 
K) and log surface gravity (from -0.5 to 5.5). This library, still 
in the visible wavelength range (3000 to 7000 A), uses the 
non-LTE model atmosphere TLUSTY (Hubenv l 19881 Hubeny 
& Lanz [T995I Lanz & Hubeny l2003t for T^s > 215QQK, 
the Kurucz's ATLAS9 model for 4750 < T^s < 21QQQK 
and Phoenix/NextGen models (Allard & Hauschildt 1995 
Hauschildt et al. .1999j which use spherical symmetry for 



cooler stars with low surface gravity. A comparison with ob- 
served spectra (from the STELIB library of Le Borgne et al. 
l2003l and the Indo-US library of Valdes et al.'2004^ shows the 
good agreement of this theoretical library with observations. 

Thus, the visible range is now quite well covered by the- 
oretical libraries, for photometric use as well as spectroscopic 
use and for a wide range of physical parameters. All the com- 
parisons with observations show that these theoretical spectra 
can reasonably well mimic real stars, at least at the spectral 
resolution where the comparisons were made. 

The goal of the present work is to go one step further, ex- 
ploring the near-infrared range where few observed fully cali- 
brated and no theoretical libraries are available. This research 
takes place in a more general framework which consists in the 
synthesis of the stellar population of galaxies hosting active 
galactic nuclei by an inverse method, described in Pelat ( I1997> 
and Moultaka & Pelat (2000 1. The H-band provides very good 
luminosity discriminators for stars later than KO (cf Dallier 
et al. 1996l and the particular region 1.57-1.64 fim of the H- 
band is clear of strong emission lines (except Brackett lines). 
It allows to sample the stellar content of the very nucleus of 
Seyfert 1 galaxies, at the contrary of the visible range where 
the strong broad emission lines of the active nucleus contami- 
nate so much the spectra of the galactic inner part that there are 
too few absorption lines from the stellar component to synthe- 
size this region. 

The lack of stellar observations at medium resolution for 
the near-infrared range, especially for super-metallic stars, 
drove us to work with theoretical spectra. But the behavior of 
model atmospheres and fluxes is not very well known in the 
infrared. 

Decin et al. (I2003> have compared several observed stars 
with theoretical spectra computed with the MARCS models 
(Gustafsson et al. [T9751 Plez et al. [T992t in the range 2.38 
to 12 yum for the ISO-SWS calibration, at a resolving power 
R ^ 1000. This study points out the difficulties of modelisa- 
tion due to strong molecular opacities and the bad accuracy and 
completeness of the atomic data in these wavelength ranges. 

In this paper, we compute theoretical spectra using the 
NeMo (Vienna New Model) grid of atmospheres (Heiter et 
al. 120021 Nendwich et al. P004 l based on the model atmo- 
sphere code ATLAS9 by Kurucz (,1993a|[T998t and Castelli et 
al. (1997 1, combined with the list of absorption lines VALD 
(Vienna Atomic Line Database, Kupka et al. T999 1, eventu- 
ally completed by molecular data collected by one of us (VT). 
These models are described in Sect.|2 Synthetic stellar spectra 
are computed with the code SynthV (built by VT), as shown 
in Sect. |3] using the model atmospheres described in Sect. |2]as 
input. Several tests on the input parameters of the spectra are 
done in Sect.|4] 

In a first step of applying our synthesis calculations 
(Sect. 13, we compare a set of observed stellar spectra with their 
corresponding models in the visible wavelength range (5000 to 
9000 A) to check the range of validity of the NeMo grid, ex- 
ploring the whole range of physical parameters (effective tem- 
perature, surface gravity and metallicity). In a second step, we 
generate synthetic spectra in the near-infrared range and com- 
pare them with observed ones. The results of this comparison 
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are described in Sect.|5] as are tests which demonstrate that the 
particular choice of model atmospheres can be expected to be 
less important than the set of line lists used for the computation 
of spectra. Our conclusions are summarized in Sect.Q 

2. Description of the model atmosplieres 

NeMo differs from the original grids of model atmospheres 
based on ATLAS9 in the treatment of the convective energy 
transport. It provides also a higher vertical resolution of the at- 
mospheres and a finer grid in effective temperature and surface 
gravity. 

This model grid of stellar atmosphere uses convection treat- 
ment without overshooting. The overshooting prescription has 
been introduced by Kurucz (I1993al I1998> and modified by 
Castelli et al. ('1997). It was supposed to take into account the 
change in the temperature gradient of the stable atmosphere 
layers near a convective zone due to the overshooting of gas 
from that zone into the stellar atmosphere. But this prescrip- 
tion is left aside in the present work, because, even if the prop- 
erties of various numerical simulations are well described and 
in good agreement when compared to observations of the Sun, 
models with overshooting are worse than models without for 
other stellar types (see Heiter et al. 120021 for a detailed discus- 
sion). 

The NeMo grids offer a choice among different convection 
models. One of them is the mixing length theory (MLT), with 
a - 0.5. The parameter a represents the ratio between the 
characteristic length (distance traveled by an element of fluid 
before its dissolution) and the scale height of the local pressure. 
This parameter is subject to discussion: according to compar- 
isons between observed and computed energy distributions for 
the Sun done by Castelli et al. J1997> . a should be set at least to 
1.25, but Van't Veer & Megessier (19961, using the same codes 
and input data as Castelli et al. ( 1997 ) . but different observa- 
tions for the Sun, found that a = 0.5 is required to fit both Hq. 
and Hjj profiles. Fuhrmann et al. (1993 ) were the first to notice 
that a value of 0.5 for the parameter a is needed to reproduce 
the Balmer line profiles of cool dwarf stars. In addition, this pa- 
rameter has to span a large domain (from 1 to 3) to reproduce 
the red giants (Stothers & Chin 1997 ). The alternative convec- 
tive models available in the NeMo grids are of "Full Spectrum 
Turbulence" (FST) type. Introduced by Canuto & Mazzitelli 
( fT99n [T992I thereafter CM model) and Canuto, Goldman & 
Mazzitelli (119961 thereafter CGM model), these models avoid 
the one-eddy approximation of MLT. In addition, both models 
were suggested to be used with a scale length different from the 
usual multiple a of the local pressure scale height (see Heiter 
et al. l2002l for further details). 

The latter models were introduced in NeMo to allow a 
choice among different treatments of the internal structure of 
the stars, depending on the aim of the model computation and 
its underlying assumption of how to describe the convective 
energy transport within the limitations of a simple convection 
model (using only algebraic rather than differential equations). 

Two levels of vertical resolution are also offered and hence 
we can either work with 72 or 288 layers. The MLT models are 



computed with 72 layers, CM models with 288 and CGM ones 
are computed for both values. 

The metallicity of the model atmopheres covers a large 
range between -2.0 and +1.0 dex and have 13 different values. 
This range of metallicity is enough for our purpose. The super 
metal rich stars, in particular, are represented with five differ- 
ent levels of metaUicity (+0.1, +0.2, +0.3, +0.5 and +1.0 dex) 
reaching the highest possible value for a real star 

NeMo provides model atmospheres for effective temper- 
atures between 4000K and lOOOOK, by successive steps of 
200K; for lowest temperatures, the model atmospheres com- 
puted with ATLAS9 become inadequate, mainly because of 
the molecular opacities which become very important for cool 
stars. The MARCS6 models (Gustafsson et al. [T975l Plez et al. 
119921 ). more dedicated to the cool stars, handle this problem 
with a more complete treatment of molecular opacity. 

The available values for the surface gravity (log g) of the 
stellar atmospheres in the NeMo grid span a range from 2.0 to 
5.0 with steps of 0.2. It is bounded at 2.0 owing to the plane- 
parallel approximation used in ATLAS9; for lower values of 
log g, spherically symmetric geometry should be used instead 
(cf. Hauschildt et al. 1999 and BaraflFe et al. 2002). 

Other models, working with the appropriate geometry like 
MARCS6 or Ph oenix/ NextGen (Allard & Hauschildt [T9951 
Hauschildt et al. I1999> are necessary for these small values. 
Indeed, MARCS6, whose main purpose is to model cool stars, 
uses also the approximation of spherically symmetric geometry 
to reproduce supergiants and the cool giants stars, which have 
a low surface gravity (Plez et al. 1992). NextGen models, like 
ATLAS9, assume LTE and plane-parallel geometry for dwarf 
stars, but a spherical symmetry is used for low-gravity giant 
and pre-main sequence stars (log g < 3.5, see Hauschildt et al. 
0999,1. 

Contrary to NeMo and MARCS6, NextGen can use a non- 
LTE model for high temperature stars. But using NLTE does 
not improve significantly the modelisation of our observed 
stars, as NLTE effects begin to occur only from 7000 K to 
higher effective temperature (Hauschildt et al. 1999), but are 
still small to at least 10000 K. Moreover, NextGen does not 
reproduce well enough the individual lines owing to the treat- 
ment of atomic and molecular lines with a direct opacity sam- 
pling method. Indeed, working with opacity distribution func- 
tions, like in ATLAS9, would ask too much computer resources 
when using NLTE calculations (Hauschildt et al. 1999). In ad- 
dition to that, too few layers are used in the published models 
to describe the bottom part of the photosphere. 

However, NextGen could be an alternative to ATLAS9 type 
model atmospheres in a next step of our project, for generating 
spectra of stars with a log g below 2.0 (for which spherical sym- 
metry is needed) and/or an effective temperature above 10000 
K. 

Bertone et al. (2004) have compared both ATLAS9 and 
NextGen models to observations in the visible range along the 
whole spectral-type sequence. The conclusions of this work are 
that both models reproduce very well the spectral energy dis- 
tribution of F type stars and earlier but this good agreement 
decreases at lower temperature, especially for K stars, owing 
to the lack of molecular treatment in those models. ATLAS9 
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provides a better fit, in general, from B to K type stars but as 
said previously NextGen is more suitable for M stars, due to the 
use of spherical geommetry for the giants and a more complete 
molecular line opacity. However, Martins et al. ( 2005 1 note that 
this comparison is made with a previous generation of NextGen 
models, using for example a mixing length parameter of 1 in- 
stead of 2, preferred by hydrodynamic models. This is also true 
for the ATLAS9 models, as Bertone et al. J2004t did not use 
the latest versions of ATLAS9, including new opacity distribu- 
tion functions (Castelli & Kurucz 2003 ), computed with more 
up-to-date solar abundances and molecular contributions than 
the previous one. 

Thus, a new comparison with observations in the visible 
range is not unnecessary. Moreover, an extensive comparison 
of spectra based on the NeMo grid of model atmospheres for 
the entire range of A to early M stars including both dwarfs and 
giants has not been done before. We hence begin our compar- 
ison with observations in the visual before proceeding to the 
infrared. The implications of changing abundances or the de- 
scription of convection at spectral resolutions relevant for stud- 
ies of galaxies are included as part of the discussion of our 
comparisons. 

3. Obtaining a theoretical spectrum 

First of all, we downloaded the model atmospheres corre- 
sponding to the stellar types wanted from the NeMo website 
(http://ams. astro. univie.ac.at/nemo/i. The models are classified 
according to the convection model (CM, CGM or MLT) and to 
the number of layers representing the atmosphere. The model 
CGM with 72 layers is detailed enough for our purpose. Models 
with 288 layers are used only for specific applications like 
the calculation of the convective scale length in stellar interior 
models (Heiter et al. l2002l l. Reduced to the medium resolution 
of our observations, both computations of a model with 72 and 
288 layers respectively give similar spectra. 

The next parameter to determine is the microturbulence ve- 
locity. For cool dwarf stars, this velocity is low: about 0-1 km/s, 
but the value is increasing towards higher luminosities, reach- 
ing values as high as 5 km/s (Gray 1992). A few stars do not 
follow this rule: hot stars, like B and O type, have a negligible 
microturbulence velocity and some specific types of A stars can 
either have a null velocity (Ap type stars, for them magnetic 
field effects are important instead) or a velocity of 4 km/s (Am 
type stars). As the microturbulence velocity has only a small 
influence on the overall shape of the spectra and on the line 
profiles at our spectral resolution, we can use a common value 
of 2 km/s for comparison with all our stellar spectra, composed 
by A to early-M type dwarf and F to K type giant stars, as 
2 km/s is a good compromise for these stars (Gray 1992 1. 

Then, the three main physical parameters of the star have 
to be chosen. The metallicity, the effective temperature and the 
surface gravity of the theoretical stellar spectrum should corre- 
spond as good as possible to the observed star to be compared. 
Therefore, once the stellar characteristics are determined, the 
nearest set of parameters (T, log g, Z) in the NeMo grid is taken. 

The metallicity is taken from Nordstroem et al. i2004j . 
Cayrel de Strobel et al. (12001^ and Barbuy & Grenon (IT990t 



when available for individual stars of the sample, otherwise as- 
sumed to be solar. 

Nordstroem et al. ('20041) have determined the effective tem- 
perature of most of the dwarf stars in our sample, for other 
stars the effective temperature is assumed, as well as the surface 
gravity, according to their spectral type from the corresponding 
values of temperature and gravity as given in Schmidt-Kaler 
(1982) and Grav( fT992l . 

Once the most suitable model atmosphere is determined, 
we can generate a theoretical flux calibrated spectrum with the 
code SynthV (by VT). This code requires several input param- 
eters such as the wavelength range for which the spectrum will 
be computed, as well as the wavelength step. This step has 
to be small enough compared to 2.5 A as the desired opacity 
in each wavelength point includes absorption from all nearest 
lines within 2.5 A, so large wavelength steps would give wrong 
results. We take 0.1 A in the visible and in the infrared range. 
Then, we can enter a rotation profile for the star, if needed, and 
indicate a list of absorption lines to be used. For our work, we 
take the Vienna Atomic Linelist Database (VALD), completed 
by several molecular line Hsts (C2, CN, CO, H2, CH, NH, OH, 
MgH, SiH, SiO, TiO, H2O '). The line profiles are approxi- 
mated by a Voigt function. SynthV also provides the possibiUty 
to change individual abundances. 

The final theoretical spectrum has to be reduced to the same 
resolution and the same sampling as the observed spectrum for 
further comparisons. So, the calculated spectrum is Gaussian 
smoothed and resampled by Fourier interpolation to the same 
step as the observed spectrum. 

4. Testing parameters 

The determination of the physical parameters of the observed 
stars is not as acurate as we would like. So it is necessary to 
investigate the nearest values of Teff / log g / Z of the grid. The 
chemical abundances can also be changed; as these abundances 
are not very well determined, it is crucial to notice how a vari- 
ation of the abundance of one element modifies the spectrum. 

The most important change is caused by the temperature. 
Indeed, the step of 200 K as was chosen for the grid computa- 
tion is still quite large for our purpose and a deviation of this 
range can be dramatic for the slope of the spectrum. The cold- 
est stars (M, K and even G type) are the ones most affected by 
a change of 200 K. 

Fig.^andlJjshow the evolution of the spectra with temper- 
ature in the visible range for the dwarfs and the giants, respec- 
tively, and Fig. |3] shows this evolution for dwarf stars in the 
infrared range. We can see that for intermediate temperature, 
there is mainly a difference in continuum. But for the extreme 
values, the modification of the spectrum is more dramatic, as it 
affects also the absorption line features. 

A change in log g causes a variation of the line profiles. 
This parameter is not very well known for observed stars. So it 
is important to test different values around a first guess. Even if 

' file VCoUjnolec.lns built by VT from Kurucz' CDROMs 15, 24, 
25, and 26, see Kurucz [l993bl and [T999l the file is available from V. 
Tsymbal upon request 
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Figure 1. Results of a variation of temperature in the visible range for dwarf stars. From the top to the bottom: T=6000K to 
4600K with a step of 200K between two spectra, arbitrarily shifted by a constant value for the purpose of clarity. 
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Figure 2. Results of a variation of temperature in the visible range for giant stars. From the top to the bottom: T=5000K to 4000K 
with a step of 200K between two spectra, arbitrarily shifted by a constant value for the purpose of clarity. 



the variation caused by a gap of 0.2 in log g is not very impor- 
tant, it can improve the comparison. 

The metallicity has an influence on the slope of the contin- 
uum: increasing the metallicity of a theoretical spectrum has a 
similar influence on the continuum as decreasing the tempera- 
ture (see e.g. Ramirez & Melendez, 2005, for a detailled dis- 
cussion). A variation of metallicity has also a clear influence 
on the strength of the absorption lines. 

The variation of individual abundances can also cause some 
changes in the spectra. SynthV uses by default solar abun- 



recent like Kurucz J1993at or Holweger J2001> give different 
values for various elements (like He, Fe, O, C, N...). These 
changes can be quite important (up to 0.2 dex for Fe). A simple 
test to probe the influence of different abundances was made 
comparing two spectra assuming the same physical parame- 
ters, except for the solar abundance of Holweger ( 2001 1 in one 
case and the one used by Kurucz Cl993 ai in the other The re- 
sult shows only slight differences, and the comparison with our 
observations cannot determine whether one is better than the 
other For our work, we chose the values given by Holweger 



dances of Anders & Grevesse il989ll. But some studies more (1200 It. 
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Figures. Results of a variation of temperature in the infrared range for dwarf stars. From the top to the bottom: T=7200K, 
T=6400K and T=6000K to 5200K with a step of 200K, arbitrarily shifted by a constant value for the purpose of clarity. 



In order to fit metallic stars, it is also important to check the 
influence of the modification of an individual element abun- 
dance on the synthetic spectrum. A change of the individual 
abundances of O and C leads to a modification of OH and 
CO line strengths, respectively. Indeed, when the ratio C/O 
increases, more CO molecules will be formed; on the other 
hand, when it decreases, more oxygen will be left to form OH 
molecules (see Decin et al.|200^. 

For the hottest stars, it is important to take into account ro- 
tational velocity. At the medium resolution of our observations, 
a convolution with a Gaussian profile is good enough to repro- 
duce this effect. Hence, we do not need to include a more accu- 
rate description of the change in the profile of the lines caused 
by the rotational velocity. Consequently, even though SynthV 
allows to compute spectra with a rotational velocity profile for 
the lines, we compute the spectra without rotational velocity in 
order to save computation time and convolve them afterwards 
with a gaussian profile. 

Additional tests have been made with a diff'erent microtur- 
bulence velocity for a very cool star (4 km/s for a MOV-type 
spectrum in the visible range) and a different number of layers 
for the convection model (288 instead of 72 for the same CGM 
model and the same physical parameters). At the resolution of 
the observed samples, these modifications do not lead to any 
difference. 



5. Results in the visible range 

Although our goal is to explore spectra in the infrared wave- 
length range, a study of the behavior of the NeMo model at- 
mospheres and spectra in the visible range provides us a good 
indication of their reliability as a function of the physical pa- 
rameters of the stars. 



5.1. Observations in the visible 

18 spectra of observed stars (A to M dwarfs and G to K giants) 
corresponding to the range of the parameters in the NeMo grid 
have been compared to theoretical stellar spectra. This sample 
of observed spectra at a resolving power of R ^ 600 is taken 
from the stellar library used by Boisson et al. <2000l l. One part 
of these observations comes from the stellar library of Silva & 
Cornell (1992i, made at the KPNO with the MARK III spec- 
trograph, they cover the wavelength range 3500-9000A. These 
spectra are, for most of them, a mean-value of several stars of 
nearby spectral type. The name of these stars and the mean as- 
sociated spectral types as given by Silva & Cornell are listed 
in Table [2 The remaining of the library, mainly supermetallic 
stars, were observed by Serote Roos et al. (Q^996 1 at the CFHT 
with the Herzberg spectrograph and at OHP with the Aurelie 
spectrograph. The spectral range is limited to 5000-9000A. The 
name, spectral type (or associated mean-spectral type) and pa- 
rameters of these stars are listed in Table[2 

The atmospheric bands are removed from the observed 
spectra to compare with the theoretical spectra. 



5.2. Comparisons 

When the physical parameters (see Table^J are known, we took 
the model with the nearest values, otherwise we used the mean 
values according to the spectral type as starting points and in- 
vestigated the nearby values to find the best agreement between 
the observed and the computed spectra. The theoretical spectra 
are computed with a wavelength step of 0. 1 A, corresponding 
to a resolution of 60000 at 6000 A, then Gaussian smoothed to 
the resolution of the observed spectra. By Fourier interpolation, 
we reduce the computed spectra to the same wavelength step as 
the observations. Spectra are normalized to 1 in the range 5440- 
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Table 1. List of observed stars in the visible, with the values of the parameters taken from the mean-values listed by Gray (I1992> 
and Schmidt-Kaler J1982> or from (1) Nordstroem et al. ( EUO?! . (2) Cayrel de Strobel et al. (2001 1 and (3) Barbuy & Grenon 
<1990> . In the last column, when no information on metallicity is available, a tick mark replaces it. The quantity < vsini > is 
given in km/s. 



Name 


Spectral Type 


< V sin I > 


reff(K) 


logg 


[Fe/H] 




HD 11 6608 ,HD 1 90785 ,HD 1 24320 


Al-3 V 


145 


8900 


4.2 


- 




& HD221741 














HD88815 


F2 V 


90 


7244' 


4.3 


-0.13' 




HD 187691, HD 149890 


F8-9 V 


7 


6026', 5902' 


4.4 


0.07- , - 


-0.44' 


HD121370 


rGOIV 


5 


5957' 


4.4 


+0.27^ 




HD38858 


G4 V 


3 


5636' 


4.5 


-0.26' 




HD161797 


rG5 IV 


3 


5700 


4.5 


+0.23" 




HD149661,HD151541,HD33278, 


G9K0 V 


2 


5176', 5236', 5300, 


4.5 


0.01', - 


-0.36',- 


HD23524,SAO66004,SAO84725 






5200', 5300, 5300 




-0.49' 




HD93800 


rKO V 


2 


5250 


4.5 


+0.43^^ 




HD39715 


rK3 V 


1 


4850 


4.6 


+0.33^^ 




HD36395 


rMl V 


1 


3850 


4.6 


+0.6^ 




HD15866, HD25894, HD2506 


GO-4 III 


10 


5500 


3.0 






HD163993 


wG8 III 


3 


4950 


2.7 


-O.P 




HD72324 


G9 III 


3 


4900 


2.7 






HD33506, HD112989 


rG9K2 III 


2 


4700 


2.6 


+0.14^ 




SAO76803 


K2 III 


2 


4500 


2.5 






HD181984 


rK2 III 


2 


4500 


2.5 


+0.P 




HD176670 


rK3 III 


1 


4300 


2.2 


-0.03^ 




HD154733, HD21110 


K4III 


1 


4000 


2.0 


-0.142 





1.2 n ' 1 ' 1 1 1 1 1 1 1 1 ' 1 ' 1 ' 1 1 r 




-0.2 LJ ' ' ' ' ' ' ' ' ' 1 ' ' ' ' 1 ' ' ' L 

5000 6000 7000 8000 

Wavelength A 



Figure 4. The observed spectrum (in black) is from a F2V type star (HD 88815), the theoretical one (in grey) is computed with 
the following parameters: T=7200K, log g=4.0, [M/H]=-0.1. Atmospheric bands are removed from the observed spectrum. The 
residual between the two spectra is (theoretical flux - observed flux)/theoretical flux. 



5460A. When the star has a rotational velocity, we convolve the 
corresponding computed spectrum with a Gaussian of the same 
velocity. 

The agreement between the observed and computed spec- 
tra is satisfactory for eff'ective temperatures ranging 4600 to 
9000 K, 9000 K corresponding to the highest temperature for 
the stars composing our sample. For these spectra, the main 



discrepancies, which consist in diff'erences in the slope of the 
blue extremity of the continuum, can be explained by the dif- 
ficulty to have a good flux calibration at the wavelength ends 
of the observational data (in particular at 5000A where strong 
Mgl, MgH and Pel absorptions are present). Various examples 
of comparisons for these stars are shown in Fig. l4l7l The spec- 
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1.2 p ' 1 ' 1 1 1 1 1 1 1 1 ' 1 ' 1 ' 1 1 r 




0.2 H 1 1 1 h- 



^ 0.1 - 




-0.2 LJ ' ' ' ' ' ' ' ' ' 1 ' ' ' ' 1 ' ' ' ^ 

5000 6000 7000 8000 

Wavelength A 

Figure 5. The observed spectrum (in black) is from a GOIV type star (HD 121370), the theoretical one (in grey) is computed with 
the following parameters: T=6000K, log g=4.4, [M/H]=+0.3. 



1.2 p < 1 < 1 1 1 1 1 1 1 1 < 1 < 1 < 1 1 r 




0.4 - 



0.2 H \ 1 \ H 



X 0.1 - 




-0.2 LJ 1 1 1 1 ^ 1 1 1 1 1 1 1 1 1 1 1 1 1 ^ 

5000 6000 7000 8000 

Wavelength A 

Figure 6. The observed spectrum (in black) is from a G8III type star (HD 163993), the theoretical one (in grey) is computed with 
the following parameters: T=5000K, log g=2.8, [M/H]=-0.1. 



tral type of the observed stars and the physical parameters used 
to compute the theoretical spectra are noticed for each figure. 

Three stars have a temperature below 4400 K, one dwarf 
and two giants. The M dwarf star (HD 36395) is particular as it 
has a very high metallicity. It is the most metallic star of Cayrel 
de Strobel's catalog.^ Indeed, to fit correctly the observations, 
we need to set the metallicity of the theoretical star as high as 
possible (Fig. 01, but a metallicity of [M/H]=+1.0 dex is not 
realistic. Fig. |5] shows that [M/H]=+0.5 dex is not sufficient. 

2 [Fe/H]=+0.6 dex 



However, recently Woolf & Wallerstein (I2005> have found a 
temperature of 3760 K (instead of 3850 K) and a metallicity 
of [Fe/H]=+0.2 dex for this star So the discrepancies may be 
simply due to the difference in temperature, out of reach for 
NeMo. 

The two others are K giant stars. First it is worth to say that 
at this low temperature, a difference of 200 K causes a vari- 
ation of the slope of the continuum more important than for 
the highest temperatures. Thus, the continuum of the K3III ob- 
served star can neither be correctly fitted by a computed spec- 
trum with a temperature of 4400 K nor with 4200 K. The first 
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Figure?. The observed spectrum (in black) is from a KOV type star (HD 93800), the theoretical one (in grey) is computed with 
the following parameters: T=5200K, log g=4.4, [M/H]=+0.5. 



K 1.6 
P 

E 




H 1 1 1 1 H 




7000 
Wavelength A 



Figure 8. The observed spectrum (in black) is from a Ml V type star (HD 36395), the theoretical one (in grey) is computed with 
the following parameters : T=4000K, log g=4.6, [M/H] =+ 1 .0. The scale for the residual is twice the scale of the previous figures. 



one is too blue and the second one is too red. Another difficulty 
comes from the Na line at 5894 A and the MgH band computed 
too strong compared to the observations. These absorption fea- 
tures are very sensitive to a variation of temperature and a too 
strong value means that the temperature is too low, but increas- 
ing the temperature would lead to a continuum too blue. This 
is particularly clear for the K4III star, which lies on the bound- 
aries of the NeMo grid for temperature and gravity (T=4000 K, 
log g=2.0), as we can see on Fig. [TO] This star comes from the 
stellar library of Silva & Cornell (1992) whose original spec- 
tra extend down to 3500A. The agreement for the slope of the 



continuum is satisfactory, but the computed Nal line, the MgH 
band as well as other lines and molecular bands such as e.g. 
Call, CN, Gband are by far too strong. If we set the tempera- 
ture at 4200 K, which is still reasonable for this kind of star, the 
accordance for the lines and molecular bands would be better, 
but the slope of the observed continuum would be too red. 

The hypothesis of plane-parallel geometry of the model be- 
gins to become unrealistic (low log gravity) and the molecu- 
lar opacities, which are not taken sufficiently into account in 
ATLAS9, become important for these cool stars. 
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H 1 1 1 1 1- 





Wavelength A 



Figure 9. Same as Fig.|8l with [M/H]=+0.5 for the theoretical spectrum. 




6000 

Wavelength A 



Figure 10. The observed spectrum (in black) is the mean-value of two K4111 type stars (HD 154733 and HD 211 10), the theoret- 
ical one (in grey) is computed with the following parameters: T=4000K, log g-2.Q, solar metallicity. The scale of this plot is not 
the same as for the previous figures. 



6. Results in the infrared range 

We can immediatly say that there are, by far, more discrepan- 
cies between the computed and the observed spectra in the in- 
frared than in the visible range, even at the medium resolution 
we have. 



6.1. Observations in the infrared 



are observations at a resolving power of R ^ 3000 at 1.6 //m 
with the KPNO Mayall 4 m Fourier Transform Spectrometer; 
these spectra have to be calibrated in flux. The stars from 
Boisson et al. l2002i) come from the ISAAC spectrograph, 
mounted on the VLT telescope, at a resolving power of R ^ 
3300 at 1.6 jum. From these samples, we selected 23 stars 
matching the available parameters of NeMo (A to M dwarfs 
and F to K giants), listed in Table|3] 



The observed spectra for this wavelength range come from 
Meyer et al. (|1998) and Boisson et al. J2002> . The Meyer's ones 
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Table 2. List of observed stars in the infrared, with the values of the parameters taken from the mean-values listed by Gray ( |1992^ 
and Schmidt-Kaler ( fTWt or from (1) Nordstroem et al. dMMl . (2) Cayrel de Strobel et al. 1200 It and (3) Barbuy & Grenon 
<1990l . In the last column, when no information on metallicity is available, a tick mark replaces it. 



Name 


Spectral Type 


< V sin i >(km/s) 


reff(K) 


logg 


[Fe/H] 


HD159217 


AO V 


150 


9700 


4.3 


- 


HD27397 


FOIV 


120 


7100 


4.3 


- 


HD48501 


F2 V 


90 


6850 


4.3 


+0.01 


HD26015 


F3 V 


65 


6776' 


4.3 


+0.11' 


HD30606 


F6 V 


10 


6152' 


4.4 


-0.01- 


HD98231 


F8.5 V 


7 


5794' 


4.4 


-0.35- 


HD112164 


rGl V 


5 


5768' 


4.4 


+0.24^ 


HD10307 


G1.5 V 


4 


5781' 


4.4 


-0.04^ 


HD98230 


G2 V 


4 


5794' 


4.5 


-0.342 


HD106116 


rG4 V 


3 


5572' 


4.5 


+0.15- 


HD20618 


G6 IV 


2 


5600 


4.5 




HD101501 


G8 V 


2 


5408' 


4.5 


+0.03- 


HD185144 


KO V 


2 


5212' 


4.5 


-0.29' 


HD22049 


K2 V 


2 


5117' 


4.6 


-0.14- 


HD39715 


rK3 V 


1 


4850 


4.6 


+0.333 


HD131977 


K4 V 


1 


4700 


4.6 


+0.03- 


HD201902 


K7 V 


1 


4100 


4.6 


-0.63- 


GL338 


MOV 


1 


3900 


4.6 




HD89025 


FO 111 


80 


7100 


3.4 




HD432 


F2 III 


75 


7278' 


3.2 


+0.18' 


HD107950 


G6 III 


4 


5050 


2.8 


-0.16- 


HD197989 


KOIII 


2 


4800 


2.7 


-0.18- 


HD3627 


K3 III 


2 


4300 


2.0 


+0.04- 










1 1 1 1 1 1 1 1 1 1 


—1 — \ — 1 — \ — 1—- 









1.57 1.58 1.59 1.6 1.61 1.62 1.63 1.64 1.65 1.66 1.67 

Wavelength 



Figure 11. The observed spectrum (in black) is from a F8.5V type star (HD 98231), the theoretical one (in grey) is computed 
with the following parameters: T=5800K, log g=4.6, [M/H]=-0.3dex. 



6.2. Comparisons 

Around 1 .6 jum hot stars are dominated by the Brackett lines, 
and a good determination of the rotational velocity of these 
stars, which broadens the lines, is very important to have the 
best possible match between the computation and the observa- 
tion. In our wavelength range, the Brackett lines at 1 .588, 1 .6 1 1 



and 1 .641 /im are nearly the only features of the observed spec- 
tra. They are well fitted by the theoretical spectra. 

When the temperature decreases, some atomic features ap- 
pear and the comparison between observed and computed spec- 
tra deteriorates. Indeed, for the F6V, a quite large amount of 
metallic lines, visible in the observed spectra, are not present or 
too weak in the computed spectra and this trend continues with 
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the F8.5V ("Fig. llll i. The continuum of these observed spectra 
is very well reproduced, but this is not the case for the lines: 
most of the metallic lines are computed too weak. 

In addition, for the F8.5V star, the Brackett lines at 1.611 
and 1.641 fim, are computed too strong. Indeed, the Brackett 
lines have almost disappeared in the observed spectrum but are 
still strong in the computation. 

These Brackett lines are also present in all the theoretical G 
stars, which is not always the case for observed stars, as seen in 
Fig- El The behaviour of the computed Brackett lines towards 
the temperature is shown on Fig.|3l the Brackett lines are still 
present in theoretical spectra for temperatures as low as 5200K. 

Then, when the temperature decreases, as expected, the 
Brackett lines are fainter and atomic lines (Fel, Sil, Mgl and 
Cal) become stronger and stronger for the observed spectra, as 
well as the theoretical spectra. But the agreement between the 
computed and the observed stars becomes worse. From Fig. 1121 
to 1151 the residuals between observed and theoretical spectra 
show that several absorption lines are missing in the theoreti- 
cal stars, iron lines for the most part. The model atmosphere is 
not the reason because the continuum shape is very good and 
several lines match perfectly, but the line list needs to be im- 
proved. 

In the infrared range, the lack of several metallic and 
molecular lines causes the discrepancies, with enhanced dif- 
ferences at low temperature due to the greatest strength of the 
lines for the coolest star (Fig.ll5>. Fi g . 1 1 31 and 1 1 41 present two 
similar stars (K2V and K4V, respectively), the first one is from 
Meyer and the second one is a VLT observation at higher res- 
olution. Both comparisons show this lack of absorption lines, 
with more details visible for Fig. 1141 

The comparison for the coolest dwarf star of this sample, a 
MOV, is not so bad for such a low temperature. The continuum 
is good in spite of the limitations of the model (Fig. I16> . We 
notice, however, that at the contrary of the previous spectra, 
the absorption lines are computed too strong for the theoretical 
spectrum, as seen thanks to the residual. This is probably due to 
the limit of validity of the model atmospheres, as already seen 
in the visible range. 

In order to investigate further which lines are missing in the 
computations, we have compared the high-resolution spectrum 
of the well-known Kllll star Arcturus (Hinkle et al. 1995) to a 
theoretical spectrum computed with the parameters of this star 
(reff = 4400 K, log g = 2.0, [M/H]=-0.2, v sin / = 3.5 km/s) and 
point out the discrepancies: Fig.^]shows a detail of this com- 
parison, and Table|3lists the missing lines in the whole range. 
In addition to the lines quoted in Table|3l several other features 
are computed too weak, in particular OH and CO molecular 
bands, certainly due to an inaccurate determination of the os- 
cillator strengths, as discussed in Lyubchik et al. J2004> . 

This study, based on the NeMo grids of atmospheres, re- 
mains valid for the entire familly of the ATLAS models. 
Indeed, as shown in Fig.EJ two theoretical spectra computed 
for the same physical parameters, with the NeMo grid and the 
ATLAS9 models with the overshooting prescription (Kurucz 
I1993al [T^98. Castelli et al. 1997 respectively), are very sim- 
ilar at our spectral resolution. The discrepancies between the 



Table 3. List of the missing lines, according to the comparison 
between Arcturus and a computed spectrum. The Ni line is not 
missing as such but shifted by 2 A in the atomic database. 



wavelength Qim) 


element 




wavelength (fim) 


element 


1.5764 


Fe 




1.6208 


Fe 


1.5893 


Fe 




1.6214 


Fe 


1.5895 


Fe 




1.6231 


Fe 


1.5913 


Fe 




1.6285 


Fe 


1.5939 


Fe 




1.6316 


Fe 


1.5954 


Fe 




1.6319 


Fe 


1.5968 


Fe 




1.6362 


Ni 


1.6007 


Fe 




1.6394 


Fe 


1.6008 


Fe 




1.6440 


Fe 


1.6041 


Fe 




1.6450 


OH 


1.6071 


Fe 




1.6517 


Fe 


1.6076 


Fe 




1.6524 


Fe 


1.6088 


Fe 




1.6532 


Fe 


1.6116 


Fe 




1.6569 


Fe 


1.6126 


Fe 








1.6175 


Fe 








1.6195 


Fe 









two different theoretical spectra are very faint compared to the 
discrepancies between the models and the observed spectra. 

The same ATLAS9 models, but without overshooting 
(NOVER models, Castelli et al. 1997), present even less dif- 
ferences with the NeMo spectra, in particular the slight dis- 
crepancy found for the Brackett lines disappear. They are more 
sensitive than other lines to the fact that the Kurucz overshoot- 
ing prescription changes the temperatures at Rosseland optical 
depths of 0.1 to 0.5. 

7. Conclusions 

In spite of some discrepancies, the comparisons between ob- 
served and theoretical spectra in the visible range suggests a 
reasonable agreement, even at the limits of the parameter range 
of NeMo. However, we have to be careful when we compute a 
model near the lower limit in temperature. 

The spectra modelized with NeMo can be used to build a 
theoretical spectral library for A to K dwarf and giant stars in 
the visible range, but this is not the case for the near-infrared 
range. 

Indeed, in the range 1.57 to 1.67 jim, the spectra com- 
puted do not reproduce very well the observations. Albeit the 
good agreement for the overall flux distribution shape, we 
can see that there are many differences for the line features 
when focusing on details of the spectra. The strength of the 
infrared absorption lines is usually underestimated in calcula- 
tions, and some lines are simply missing (Fe, OH and CO lines 
are the most problematic ones). As pointed out by Decin et al. 
12003) for the MARCS6 models and Lyubchik et al. 
for NextGen models of ultracool dwarfs and a Kurucz model 
for Arcturus, it was not possible to generate synthetic spectra 
which can reproduce observed spectra in the infrared with the 
line lists that have been used in constructing the model atmo- 
spheres, even at a medium resolution. In particular, the oscilla- 
tor strengths are still not known sufficiently well. 
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1.57 1.58 1.59 1 



Figure 12. The observed spectrum (in black) is from a G4V type star (HD 106116), the theoretical one (in grey) is computed 
with the following parameters; T=5800K, log g=4.4, [M/H]=+0.3dex. 




H — ' — h 



H — ' — \ — H 



H \ ' h 



H — H 




1.57 1.58 1.59 1.6 1.61 1.62 1.63 1.64 1.65 1.66 1.67 

Wavelength fj.m 

Figure 13. The observed spectrum (in black) is from a K2V type star (HD 22049), the theoretical one (in grey) is computed with 
the following parameters: T=4800K, log g=4.6, [M/H]=-0.1dex. 



We have also performed a comparison of spectra that in- 
clude the lines from iron peak elements with predicted energy 
levels, as published by Kurucz (1998), with the observations 
of Arcturus discussed above. As to be expected, such spectra 
contain more lines, but the inaccuracy of their energy levels 
frequently places them at the wrong wavelengths and the over- 
all match hardly improves (the total flux distribution over all 
wavelength ranges, particularly in the ultraviolet, is closer to 
observations when including this set of lines, but for the lim- 
ited wavelength range around the H band the effects are small 
and sufficiently compensated when setting the zero point of the 



flux distribution). For the case of Arcturus we also performed 
a comparison with spectra computed with PHOENIX in LTE 
(P. Hauschildt, priv. comm. 2005). The resulting spectra for the 
1.57 to 1.67 jum range were found to be rather similar to those 
from NeMo/VALD/SynthV when including the predicted level 
lines. One important reason for this is certainly the fact that the 
atomic line lists for PHOENIX are essentially those of Kurucz 
|1998). Because the flux distribution of the PHOENIX spectra 
is similar to the observations as well, at least for the K giants 
the detailed choice of the model atmosphere code appears to be 
clearly less important than the choice of atomic line Usts (note 
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1.57 1.58 1.59 1 



Figure 14. The observed spectrum (in black) is from a K4V type star (HD 131977), the theoretical one (in grey) is computed 
with the following parameters; T=4800K, log g=4.6, solar metaUicity. 





1.57 1.58 1.59 l.e 



1.61 1.62 1.63 1.64 1.65 1.66 1.67 

Wavelength f^m 



Figure 15. The observed spectrum (in black) is from a K3111 type star (HD 3627), the theoretical one (in grey) is computed with 
the following parameters: T=4400K, log g=2.0, solar metallicity. 



that PHOENIX uses its own collection of molecular line lists, 
different from the one we have used here). Considering the uni- 
formity of the deterioration of the match of spectra in the 1.57 
to 1.67 jum range when looking at the sequence from F to K 
stars we conclude that the insufficient line lists, and in particu- 
lar lists of atomic lines, are the main obstacle for a more satis- 
factory match of observed spectra of these groups of stars. The 
modelizing in the infrared range needs some further improve- 
ments, in particular for the absorption lines database, before to 
build a theoretical spectral library which can be used with high 
benefit instead of an observed star library. 



The lack of M stars in spectral library would be very much 
prejudicial to the study of stellar populations as the variations 
of their strong atomic lines and molecular bands along their 
evolution from dwarf to supergiant to giant provides very good 
age discriminators (from 10^ to lO'^yrs). M stars peak in a 
wavelength range which is not much absorbed even in heavily 
reddened region, as young stellar clusters, making them easily 
detectable. Moreover they are known to be very important con- 
tributors to the stellar populations of galaxies as well for the 
mass as for the luminosity following the age of the population. 
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Figure 16. The observed spectrum (in black) is from a MOV type star (GL 338), the theoretical one (in grey) is computed with 
the following parameters: T=4000K, log g=4.6, solar metallicity. 




1.57 1.58 1.59 1.6 1.61 1.62 1.63 1.64 1.65 1.66 1.67 

Wavelength fj.m 

Figure 17. Comparison between two theoretical models (ATLAS9 and NeMo) with the following parameters: T=4400K, log 
g=2.0, solar metallicity. Note that the scale for the residual is enhanced compared to all other figures. 



All this make a good theoretical library of M stars very critical 
in order to extend incomplete observed library. 

The prospects of constructing such a library from 
the upcoming generation of model atmospheres (MARCS, 
PHOENIX, perhaps future versions of ATLAS) are indeed im- 
proving because of the enormous efforts spent in extending the 
molecular Une data and equation of state. To match the spectra 
of the hot end of M stars in the H band will nevertheless require 
more complete atomic data, although this is less crucial as for 
K stars. Efforts along this direction are cmrently made. 
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1.57 1.58 1.59 1.6 1.61 1.62 1.63 1.64 1.65 1.66 1.67 



Figure 18. Comparison between a high resolution observation ^ 100000) of Arcturus (in black) and the corresponding 
theoretical star (in grey). The flux is given with continuum normalized to 1 in order to better show the missing lines. 
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Figure 19. The same as Fig.[^but zoomed in into a Umited wavelength range. 



the importance of the line lists used relative to the particular choice of 
model atmosphere codes. 
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